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PRODUCT REVIEW 


printed, faxed, or saved to a Post- 
Script file. Not only has printing 
become a nonissue in our NeXTSTEP 
programs, but adding windows to an 
application for expository text or 
pictures is trivial. Drag a window 
from the palette in IB, drag a text 
object into the window, and start 
typing. The window will even create 
scroll bars if the text exceeds the size 
of the exposed view. 


Serious contender 


With the exception of a good Pascal 
environment for teaching our intro- 
ductory computational-physics 
course, we have been able to fill every 
departmental application need either 


through public domain, shareware, or 
commercial applications at PC prices. 
Our current standard application 
suite includes word processors, sym- 
bolic-computation package, equation 
editor, image editor and drawing 
tools, and a spreadsheet. 

NeXT ought to be considered a 
serious operating system for compu- 
tational physics, particularly in small 
departments with limited computer 
support. NeXT bundled applications 
fill many system-administration, 
communication, and programming 
needs. Reasonably priced commercial 
applications exist. Once NeXTSTEP 
has been set up and is running, it 
requires less software administrative 


work, and is far more functional, than 
single-user operating systems. Our 
students have given NeXTSTEP the 
highest compliment—when identical 
applications are available on two or 
three different operating systems, 
they usually choose to work with 
NeXTSTEP. 

This review is based in the 
author’s experience managing and 
programming a cluster of 14 NeXT- 
stations, as part of a NeXT, DOS- 
/Windows, Apple network in the 
Davidson College Physics Depart- 
ment. This department’s NeXTSTEP 
program LASER won an honorable 
mention in the 1992 Computers in 
Physics software competition. 


INTERACTIVE PHYSICS II: A PHYSICS 
SIMULATION LABORATORY FOR THE MACINTOSH 


A. John Mallinckrodt 


interactive Physics Il system re- 
quirements: 

Macintosh computer with System 
6.0.5 or later, hard disk, 1500- 
kbytes free RAM for black and 
white or 2500-kbytes free RAM for 
8-bit color. 

Available from: Knowledge Revolu- 
tion, 15 Brush Place, San Francis- 
co, CA 94103. Order no. (800) 766- 
6615; Tel: (415) 553-8153; Fax: 
(415) 553-8012. List price: IP Il 
$399; IP Il student version (no 
technical support or upgrade rights) 
$99; Fun Physics (the original Inter- 
active Physics) $99; educational 
discounts available for lab packs. 


John Mallinckrodt is Associate Professor of 
Physics at California State Polytechnic Uni- 
versity, Pomona, CA 91768. He likes to call 
himself a generalist; he is, however, particu- 
larly interested in the use of computer-based 
simulations in physics education. E-mail: 
ajmallinckro @ csupomona.edu. 


N2 quite four years ago, at the 
summer meeting of the Ameri- 


can Association of Physics Teachers 
in San Luis Obispo, CA, attendees 
were buzzing about a remarkable new 
Macintosh application being demon- 
strated by its young author, David 
Baszucki. By now virtually every 
physicist with an interest in educa- 
tional software is aware of that prod- 
uct—Interactive Physics (IP)—and 
its unique ability to simulate the 
dynamics of two-dimensional me- 
chanical systems that are literally 
drawn using the mouse. The original 
package, reviewed in the May/June 
1991 issue of Computers in Physics," 
has probably had more impact upon 
the use of computers in primary 
physics instruction than any other 
single application, with the possible 
exception of spreadsheets. Last sum- 
mer Knowledge Revolution unveiled 
the greatly enhanced Interactive 
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Physics II (IP II). In this review I 
shall discuss only a few of the more 
significant and advanced features 
making their debut in IP II. 

While IP II looks and works 
very much like its predecessor, and 
runs experiments written for IP, the 
underlying “engine” and data struc- 
tures have been rewritten in order to 
build a more flexible and reliable base 
for future improvements. IP II also 
introduces an impressive array of new 
features: 


e User-definable forces that take 
the form of two-body or field interac- 
tions. (IP II includes predefined 
forces simulating planetary and ter- 
restrial gravity, electrostatic and 
magnetic fields, and air resistance.) 


e Customizable integration al- 
gorithms (see below). 


e Actuators and motors that can 
exert user-defined forces or enforce 


user-defined relative positions, veloc- 
ities, and accelerations. 


@A large number of new con- 
straints including rotational springs 
and dashpots, rods, separators, piv- 
ots, keyed or unkeyed slot joints, and 
pulley systems. (Note: Pulley systems 
consist of ropes constrained to pass 
without friction through one or more 
fixed points.) 


e “Controls” that allow the user 
to change parameters of an experi- 
ment during a run and without the 
need to invoke dialog boxes. 


@ User-definable graphs and me- 
ters. 


e Greatly enhanced control over 
the properties and appearances of 
objects, including the ability to im- 
port pictures (from painting and 
drawing programs) and, optionally, 
to attach them to objects (see Fig. 1). 


e@An improved algorithm for 
handling frictional contact between 
objects, which models both static and 
kinetic friction. 

eThe ability to view experi- 
ments from reference frames—iner- 
tial or otherwise—attached to any 
object or the system center of mass. 


è “Player documents” that 
shield users from the program’s full 
complexity and give them access to a 
limited number of experimental pa- 
rameters via “menu buttons” and 
“controls” (see below). 


@ Quicktime output of experi- 
mental runs for more rapid playback. 


But by far the most important 
addition, and the one that underlies 
and enhances the flexibility of many 
of those already mentioned, is a new 
formula language that allows users to 
access simulation parameters and to 
use the values of those parameters to 
control others. 


Formulas in IP Il 


The properties and appearance of 
objects in IP II are specified in the 
text fields of a pair of floating win- 
dows; for example, one can specify 
the initial velocity of a body, the 
torsion constant of a rotational 
spring, or the variable(s) to be plot- 
ted in a graph. Much of the power of 
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Fig. 1: Frame from an IP II simulation shows the use of imported graphics, controls, menu but- 


tons, formatted text, and meters. 


IP II, however, derives from the 
user’s ability to substitute formulas 
for these values in a manner analo- 
gous to the use of formulas in spread- 
sheets. 

The IP II formula language 
assigns to every object an identifier, 
and to most identifiers a variety of 
fields that contain information about 
the object’s properties and current 
dynamic state. Through these the 
user has access to the values of most 
simulation parameters, including 
constraint forces and lengths, meter 
values, control outputs, friction and 
normal forces, and the current frame 
number. One combines or operates on 
these values with logical and/or 
mathematical functions to form for- 
mulas that control the properties of 
other objects. Thus, it is a relatively 
simple matter to create a meter that 
measures, for instance, the angular 
momentum of a body relative to some 
point, moving or not, simply by 
inserting the proper formula into the 
field that describes the output of a 
meter. In a more advanced applica- 
tion, one can control the strength of a 
magnetic field according to what time 
it is, according to the location of an 
object, or even, if one is so perversely 


inclined, according to the tension in a 
pulley system. Indeed, if one cares 
to—or is simply careless—one can 
quite easily undermine the prepro- 
grammed intention of the IP II 
engine, which is, of course, to simu- 
late a world governed by Newton’s 
laws. 

There are a few inconsistencies 
in the syntax of the formula language, 
some confusion and errors in the 
manual descriptions of functions, and 
a few mistakes in the formulas that 
describe the default measurements of 
some meters. For example, some 
forces treat the x and y components of 
a vector as being along and perpen- 
dicular to the line between two ob- 
jects instead of horizontal and verti- 
cal. This is often desirable, but it is 
not always obvious which rule ap- 
plies, and the manual even gets it 
wrong in at least one case—that of 
two-body constraint forces. I suspect 
that it is at least partially because of 
this confusion that the default formu- 
las in meters measuring the gravita- 
tional or electrostatic potential ener- 
gy of two bodies are wrong and, as a 
result, give values that depend upon 
the direction from one body to the 
other. 
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PRODUCT REVIEW 


Integration algorithms 


At the heart of IP II is the “engine” 
that carries out the‘discrete integra- 
tion of Newton’s second law. IP II 
allows the user to modify the integra- 
tion algorithm in several ways: One 
has the choice of three basic meth- 
ods—Euler, fourth-order Runge- 
Kutta, and predictor-corrector. One 
can choose fixed or variable simula- 
tion time steps and can also specify a 
longer time step between screen up- 
dates than is used for the simulation 
itself. (This is a handy feature for 
increasing the accuracy of the simula- 
tion without unnecessarily increasing 
the graphics burden on the proces- 
sor.) One can ask the engine to 
produce warnings when it detects 
high velocities and/or accelerations 
or constraints that are redundant or 
inconsistent. Finally, there is an im- 
portant but not too well-explained 
option that pertains to a fundamental 
difficulty any discrete integrator will 
have dealing with collisions, especial- 
ly when those collisions involve mul- 
tiple bodies, frictional contact, and 
partial elasticity. Because the han- 
dling of such collisions involves a 
number of somewhat ad hoc opera- 
tions, neither option here is wholly 
satisfactory. In general, the user is 


well advised to view the detailed 
results of collision intensive experi- 
ments with appropriate skepticism. 

I tested the integration algor- 
ithms with the following experiment: 
A 1.0-kg body is attached to a damper 
that applies a force given by —kv° with 
k = 1 N/(m/s)?. The body is given 
an initial velocity of 10 m/s. This 
problem is somewhat challenging to 
simulate because the retarding force 
is initially so large (1000 N) and 
decreases so quickly. Note, for in- 
stance, that using Euler’s method 
with a fixed time step of more than 
0.01 s will yield a retarding impulse in 


the first time step that is greater than ` 


the initial momentum of the body. As 
a result, the body’s direction of 
motion will be immediately re- 
versed—an unphysical result of the 
overly large time step. Using a 
smaller time step would seem to be 
necessary to obtain reasonable re- 
sults. On the other hand, as the 
velocity decreases, the per-step im- 
pulse falls rapidly and a smaller time 
step is computationally inefficient. 
Problems like this benefit greatly 
from the use of variable-time-step 
methods. 

The table shows the results of the 
six integration methods provided by 


IP II applied to this problem. For a 
range of basic step sizes, the table 
shows the time (in seconds) that it 
takes the body to move 2.0 m, as 
measured both in simulation time 
(for which an analytical calculation 
yields the result 2.20 s) and in real 
time on a Mac IIci running in 16- 
color mode. As you can see, for this 
problem, the variable-time-step pre- 
dictor-corrector method was the 
stand-out champion, giving accurate 
results for any time step up to the 
maximum used. As expected, the 
fixed-time-step Euler method fails for 
time steps of 0.01 s or more. 

As another test of the algorithms 
I created a simulation of a satellite 
orbiting a planet in a highly eccentric 
orbit. This kind of simulation often 
causes problems because of the large 
spatial variations in the force and the 
high velocities that occur near the 
point of closest approach. I measured 
the orbital period, the percentage 
change in energy during the first 
orbit, and the computation time re- 
quired to complete the orbit. In this 
case the variable-time-step Runge- 
Kutta method emerged as the best 
choice, with the predictor-corrector 
method the big loser. Interestingly, 
the Euler method tended to give the 


Table. A comparison of simulation times and real times observed using a test experiment as described in the text. 
Each row gives the results of a different integration method (“PC”=predictor-corrector, “RK” =fourth-order 
Runge-Kutta, “var” —variable time step, “fix” —fixed time step), and the columns indicate the basic time steps 
used. Each entry gives “simulation time/real time” (both in seconds) on a Mac lici. 


0.2 s 0.1 s 


0.05 s 0.02 s 0.01 s 


0.005 s 0.002 s 


Euler var 2.80/3 2.60/5 2.50/7 2.40/15 2.360/26 2.300/51 2.238/122 


Euler fix 2.435/50 2.276/108 


PC var 2.20/3* 2.20/5 2.20/7 2.20/15 2.19/29 2.195/57 2.200/140 


PC fix 0.64/6 2.24/24 2.185/44 2.198/111 


RK var 2.20/20 2.21/37 2.205/73 2.202/182 


RK fix 2.48/16 2.26/27 2.205/52 2.202/142 


* Had to ignore a high-force warning. 
— Simulation failed with these parameters. 
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Fig. 2: Simulation of the inverted pendulum, from which this frame was taken, exercises many 
of the tools that have been introduced in IP Il. 


smallest cumulative energy change, 
even though it did vary during the 
orbit by as much as 60%! 


An example 


The inverted pendulum’ is a nonlin- 
ear system that exhibits both counter- 
intuitive dynamic stability and chao- 
tic dynamics. Fig. 2 shows a frame 
from my IP II simulation of the 
system, which uses many of the new 
tools in IP II. The rod is connected by 
a pivot at its lower end to an invisible 
object whose y coordinate has been 
specified (using the formula tan- 
guage) as a sinusoidal function of ` 
time; the formula references the out- 
put of the “Amplitude” control. The 
initial x and y positions of the circular 
body attached to the upper end of the 
rod are determined by formulas that 
reference the “Initial Angle” control. 
The circle is also subject to gravity 
and a linear velocity-dependent drag 
force that are both implemented using 
a custom force field. The formulas 
that define that field reference the 
outputs of the other two controls. 
(Note: In order to preserve the 
accuracy of the chosen time step, the 
“Frequency” of the oscillator was 
used to renormalize the gravitation 
and drag forces.) To create the phase 
plot, references to the x components 
of the circle’s position and velocity 
were inserted into the fields that 
specify the x and y coordinates of the 
graph. The simulation took about half 
an hour to create and debug, and 
reproduces many specific behaviors 


discussed in the previously referenced 
article. 


Player documents and course 
ancillaries 


IP II has two modes of operation. 
When launched it enters “Edit 
Mode” by default, giving the user full 
access to all editing tools and menu 
items. However, experiments can be 
created and then saved in “Player 
Mode.” Player mode documents 
shield the user from all commands 
except those that run, stop, and reset 
the experiment, and any others that 
the creator chooses to allow in the 
form of controls and/or menu but- 
tons. This facility encourages the 
creation of exercises that focus the 
attention of the student on the effects 
of varying specific parameters in a 
given situation. Knowledge Revolu- 
tion is promoting the use of these 
Player Mode documents in hopes of 
cultivating a lively aftermarket in the 
exchange of high-quality simulations 
produced by faculty for use in the 
classroom and for class assignments. 
Saunders College Publishing is using 
Player Mode documents to support 
its offerings in introductory physics. 

In a related vein, Knowledge 
Revolution has produced the Interac- 
tive Physics Player, and has given 
Prentice-Hall exclusive rights to its 
use. Course ancillaries, in the form of 
the IP Player and Player files keyed to 
specific examples and problems, are 
available to students using Prentice- 
Hall’s introductory physics texts. It 
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should be noted that the IP Player is 
designed to work only with the 
simulations that come with the Pren- 
tice-Hall textbooks, and not with files 
created by IP II. 


Final words 


With all of its new features, IP II is 
slower than its predecessor; on identi- 
cal simulations, I found that it typi- 
cally ran at about half the speed of IP. 
On the other hand, a judicious choice 
of integration algorithm and step size 
can counter, and even reverse, this 
effect without sacrificing accuracy. 
Furthermore, the new version ex- 
tracts smaller relative penalties for 
running in 16-color mode rather than 
black and white. And finally, it is 
robust; despite some pretty severe 
torture at my hands, the program 


features and improved accuracy. 
There are also plans for a Windows 
version. 

This review is the product of 
several months of fairly detailed prob- 
ing into the behavior of IP II and 
numerous discussions of my findings 
with David and Greg Baszucki at 
Knowledge Revolution. Early in my 
investigations, Knowledge Revolu- 
tion released a maintenance upgrade 
for IP II that addressed several of the 
problems that I and others had 
discovered. The free distribution of 
this upgrade to registered users of IP 
II demonstrates the company’s com- 
mitment to eliminating the inevitable 
bugs that crop up in a major revision. 
That job, however, is not over. There 
are still a number of annoyances, 
peculiarities, and bugs—mostly mi- 


nor, but a few more substantial. They 
tend to be hidden in the more ad- 
vanced areas of the program, where 
many users may never encounter 
them. I have sent my own detailed list 
to Knowledge Revolution and, if the 
past is any indication, most will be 
addressed in a future maintenance 
upgrade. 

Let there be no doubt about my 
affection for IP II. It is a flexible and 
fascinating product that offers a sig- 
nificantly expanded array of features 
over its predecessor. Despite its occa- 
sional idiosyncrasies, I find it to be an 
efficient tool for preparing and pre- 
senting simulations that are effective 
and readily modified in the class- 
room; for testing, sharpening, and 
often correcting one’s intuition; and 
simply for casual explorations. 


locked up only once. 

It is worth noting here that, by 
the time this review is published, 
Knowledge Revolution intends to 
have released a “professional” ver- 
sion of IP that makes use of the math 
coprocessor, and sports additional 
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